Stall, spiculate or runaway - the fate of fibers growing towards fluctuating membranes 
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We solve the dynamic equations of motion for a growing semi-flexible polymer, or fiber, approaching a 
fluctuating membrane at an angle. At late times we find three different regimes: fiber stalling, when fiber 
growth stops due to membrane resistence, run-away, in which the polymer bends away from the membrane, 
and another regime in which the membrane response is nonlinear and tubular membrane spicules are formed. 
We discuss which regions of the resulting 'phase diagram' are explored by (i) single and bundled actin fibers 
in living cells, (ii) sickle hemoglobin fibers in red blood cells, and (iii) microtubules growing within artificial 
vesicles. We complement our analysis with full 3 -dimensional stochastic simulations. 
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Semi-flexible polymer fibers, making up the cytoskeleton, 
interact with fluctating membranes continuously in living cells 
0] . The polymerisation of filaments such as actin, with persis- 
tence lengths l p ~ 10/im, can deform the outer plasma mem- 
brane of the cell, giving rise to protrusions that can be sheet- 
like (lamellipodia) or more localised (filopodia) 0]. Forces 
are thought to arise through polymerisation at the fiber tips 
Q, 01 . A similar mechanism is thought to drive growth of 
sickle hemoglobin fibers, leading to a pathological rigidifica- 
tion of red blood cells 0], as well as microtubules, such as 
can be grown in artificial vesicles 0] . 

While the phenomenon of Euler buckling of a fiber of fixed 
length impinging on a solid obstacle at an angle is well studied 
(see e.g. and Refs. therein), the case in which the elastic 
fiber is growing, and the obstacle is a fluctuating membrane, 
have so far received much less attention, despite its more di- 
rect relevance to biological systems. Thus here we propose a 
set of equations of motion that describe the coupled dynam- 
ics of a growing semi-flexible polymer close to a fluctuating 
membrane (Fig 1 ). We complement this treatment with full 3- 
dimensional dynamic Monte-Carlo (3dMC) simulations. Our 
main focus is on the dynamic regimes attained at late times. 
We find that fibers may (i) stall, (ii) run-away, i.e. bend away 
from the membrane, or (iii) cause the formation of tubular 
membrane spicules. What differentiates our approach from 
earlier work is the explicit description of the flexibility of both 
the membrane and the polymer QSSHS. which allows us 
to include all three fates of the fiber. 

The Hamiltonians governing the membrane and the poly- 
mer elasticity, respectively H m and H p , are given by : 



H m = \J d 2 r(n m (y\u) 2 ^cj(y^u) 2 ) 
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where K m and a denote the membrane rigidity and surface 
tension respectively and are typically n m ~ 10 _19 J and 
a ~ 10 -4 Jm -2 for biomembranes, u measures the normal 
deviation of the membrane from local flatness, which we take 
to be the plane z = 0, k = ksTlp = l p //3 is the polymer 



bending rigidity, L is its (instantaneous) length and 6(s,t) is 
the angle between the local direction of the fiber at arc length 
position s and the plane z = at time t, see Fig 1. Our use 
of an explicit Hamiltonian to describe the flexibility of both 
the membrane and polymer differentiates our approach from 
earlier work 0, IS 01 • The probability distribution of 
the membrane displacement can be calculated from the con- 
ditional partition function that it is at fixed height z m at one 
point. This is here approximated by JTcUl Z m = exp(— Az^), 
where A = — 7 27rcr v , ft being the area of the "frame" 

log (1+^2) 

supporting the membrane patch under consideration. A mea- 
sures the membrane 'softness' (A — > 00 corresponds to a hard 
wall), and typically A ~ 10 -4 Jm -2 fioll . The partition func- 
tion Z m can be integrated over all allowed z m to obtain the 
normal force acting on the fiber tip due simply to the steric 
repulsion between tip and membrane. As is common in stud- 
ies of buckling we perform a saddle point approximation on 
H p in Eq. [I] and consider the adiabatic limit in which, at 
all times, the fiber is in equilibrium. This is reasonable in 
all physical cases as the fiber relaxation time (at most ms) is 
much smaller than the inverse polymerisation rate (typically 
0.01-1 s). We seek to solve the following coupled equations 
for the fiber shape, growth rate and tip height and the normal 
force on the tip, respectively 
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Here /c on ,off are, respectively, the rates of polymerisation and 
depolymerisation at the fiber tip, d and A the initial and in- 
stantaneous distances between the fiber tip and the membrane 
frame, 5 the increase in length upon addition of a monomer, 
while /(A) is the force exerted by the membrane on the fiber 
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FIG. 1 : Sketch of the system under consideration. A growing fiber is 
obliquely incident onto a fluctuating membrane, which is anchored 
to some distant frame. The fiber is clamped at s — and its length 
at time t is L(t). 
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FIG. 2: Plot of the critical angle C (solid line) below which a run- 
away transition occurs, as a function of k n/k f£, for actin-like pa- 
rameters (d 100 nm, k = 10 k B T /im, VA5 - 0.8, 8 = 2.5 
nm). (The run-away state shown on the right is obtained with 3dMC 
simulations with 0o = 45°, k on /k Q E — oo and VA5 = 0.4.) 

tip fioll . We assume constant rates /c on ,off but time varying 
rates could be included to model, e.g. actin concentration 
fluctuations. The fiber is treated as if clamped at one end 
0(s = 0,t) = 0(s = 0) = 0o, into the cytoskeletal mesh, 
while the end at s = L(t) is free, (d0(s,t)/ds) 3=L , t \ = 0. 
Eq. |3 has been derived by assuming that the membrane and 
the fiber only interact via excluded volume and that this can 
be introduced at a single point (the fiber tip) fioll . This may 
fail for a highly bent fiber that often contacts the membrane 
far from its tip full . We solve Eqs. 12141 numerically, via a 
standard Euler relaxation algorithm and finite difference dis- 
cretisation. The solutions for the system of Eqs. 1215 1 fall into 
two classes. The fiber can either (a) grow until it stalls, or (b) 
undergo a run-away transition 01. This is distinct from usual 
fiber buckling, which is a purely equilibrium phenomenon. 
The physical parameters determining the system dynamics are 

the four dimensionless quantities yr^r^ ^on/^off, ksTd 2 9 
and d/S. Eq. |2| neglects fluctuations around the average 
fiber shape (see below for a discussion of the effect of these). 
Throughout we work at physiological temperature T. 

Fig. 2 shows a cut of the 'phase diagram': for 0o > C 
the fiber ultimately stalls, otherwise it is bent away by the 
membrane. Here C is established by variation of k on /k Q ^, 
for parameters typical of actin fibers in cells (see caption). 

As the incidence angle approaches the critical threshold 
from below, the growth of the fiber is strongly non-linear at 
intermediate times (see Fig. 3a). This is because the growth 
rate is reduced by a factor exponential in /(A) (see Eq. (2)). 
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FIG. 3: For three different initial incidence angle 0o we show (a) 
the fiber length as a function of time, and (b) the tangential adiabatic 
force ft as a function of length. Here k n/k E = 5 (other param- 
eters as Fig. 2), so that the stalling force is log (5)/cbT/£, and the 
fiber is stalled for 0o = 75° and runs away at the two smaller inci- 
dence angles. This diagram serves to show that all stalled states are 
metastable. The analytical prediction for the adiabatic longitudinal 
force (dotted line) is compared to the numerical result for 0o = 45°. 
The time needed to tunnel the barrier to run-away growth can be 
estimated as being proportional to the exponential of an energy bar- 
rier, being the area of the coloured region for 0o = 75° in (b). The 
linear-logarithmic plot of this tunneling time as a function of its ini- 
tial incidence angle, found by an application of Kramers' theory, is 
shown in (c). From top to bottom, k on /k f? = 1.1, 1.5, 3, 5 and 10. 



This force bends the tip (Eq. (1)) and as a result the growth 
speeds up again, as the factor sin(0(L(£), t)) decreases (Eq. 
(2)). Far from the transition line this behaviour is not found 
and the fiber growth is nearly linear at all times, see Fig. 3a. 

The physics of the transition between the metastable stalled 
state and the run-away state is controlled by the tangential 
force on the fiber tip ft(L), see Fig. 3b - f t is related to 
the normal force in Eq. |5]via a factor sin(0(L)) to be de- 
termined self consistently from Eqs. 12151 If the thermody- 
namic stall force f sta ii = (k B T/S)\og(k on /k of f) is above 
the maximum of f t the fiber runs away, otherwise the fiber 
stalls when it reaches the length at which the membrane gen- 
erated force equals its thermodynamic stalling force. How- 
ever, as the adiabatic force vanishes for very long fibers, the 
stalled state is, in fact, metastable and the fiber may still run 
away provided it can 'tunnel' through the free energy bar- 
rier present. Using elementary transition state theory, we 
can estimate the time required for such a tunneling event as 

At = To exp ^ J^ 2 dL f^j where To ~ ms is a typical mi- 
croscopic relaxation time (or inverse attempt frequency) and 
L\ and L2 are given by intercepts of the force curve with the 
thermodynamic stall force, see Fig. 3b. Fig. 3c shows a rep- 
resentative plot of At versus 0o, with actin-like parameters. 

Actin fibers in cells can be biochemically capped Q] . If the 
tunneling time is larger than the inverse capping rate, typically 
< 1 s, the fiber growth will be halted. Tunneling from the 
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metastable stalled state may thus often be negligible . 

While Eqs. I2l5l have to be solved numerically, one can em- 
ploy a Gaussian approximation for the fiber tip fluctuations, 
controlled by H p in {0, which permits some analytic analysis 
and, in particular, yields an explicit formula for ft(L). 

To achieve this, starting from Eq. [l] we adopt a standard 
path integral method IKJ1 . and compute the first two moments 
of the probability distribution for the fiber tip position 

(z p ) = 2ftsin<9 (l - exp ( - ^-)) (6) 
/ 2\ , 2 ( L ( L \ , cos 2^o ,~ 

(4)o = 4 *% + -p(-^)-i- — ^ 

(3 + exp(-f)-4exp(-|))) 

where the subscript denotes the fact that these are derived 
in the absence of a membrane. Using these moments we 
can approximate the distribution of heights z p of the tip of 
the semiflexible polymer by a Gaussian distribution, Z p = 

exp ( - B (z p - (z p ) ) 2 ) where B = \ —L^ . This 
approach is appropriate provided the force on the fiber is <C 
the Euler buckling threshold, 7r 2 ^/(4L 2 ). So far we have cal- 
culated separately : (i) the membrane deformation, and (ii) the 
fiber tip distributions. We now introduce the steric interaction 
between the rod and membrane as follows : 

/co poo 
dz m / dz p lAi^Zm — Zp) Z rn Z p 
-co J — CO 

= Vl-erff } Z ^~ d )) (8) 

where U{x) = 1 if x > and otherwise. 

It is reassuring that for completely rigid fibers we recover 
the corresponding result derived in Ref. fioll . Using Z tot 
one can calculate, e.g. the tangential force acting on the tip 
jh = — dlog ^ tot " ) . This can be directly compared with the 
full numerical treatment (see Fig. 3b, in which the compari- 
son is shown for a fiber with Oo = 45°). For small incidence 
angle, up to Oo ~ 55° (with the parameter values of Figs. 2-3), 
the agreement is good, breaking down gradually as the initial 
incidence angle approaches normal incidence. The smoother 
force rise in the analytics for = 55° is due to the inclu- 
sion of tip fluctuations, which do not alter the value of the 
maximum. (We note that if the membrane is floppier than the 
agreement persists to a larger angle.) To see why the analytics 
eventually break down, we note that invoking a Gaussian tip 
distribution amounts to replacing the fiber with a spring with 
a Hookean constant B, which becomes stiff er with 6 . For 
not too close to 90°, B ~ 2L3 ^tjqs , in agreement with the 
treatment of Ref. |6]. Our combined treatment of fiber tip 
and membrane fluctuations using a Gaussian distribution can 
be seen as the natural generalisation of the elastic fiber model 
of Refs. yllal to the case of a fluctuating soft membrane. As 
the incidence angle becomes close to normal, intuition sug- 
gests that the tip fluctuations will become 'one-sided' and the 
Gaussian approximation may then become poor. 



We note from Fig. 2 that in an eukaryotic cell (where 
&on Aoff ~ 100 01), C ~ 86° for an actin fiber. Thus the 
majority of fibers in a cytoskeletal actin network may be ex- 
pected to eventually run away from the membrane. Fibers stall 
at smaller angles than this if the initial fiber-membrane separa- 
tion is smaller (than d = 100 nm). Electron microscopy of the 
actin network 11211 reveals that some of the fibers do appear to 
be bent into a sub-membrane 'thatch work', which may help 
the cytoskeleton sustain stress although might contribute less 
efficiently to motility. However in lamellipodia and filopo- 
dia, observed at the leading edge of a moving cell, the fibres 
are actively pushing with their tips pointing towards the mem- 
brane 1 13]. It is then natural to ask what 'countermeasures' a 
moving cell takes to avoid run away. 

One efficient strategy is to combine many fibers into a 
thicker, thus stiffer, bundle. If n fibers are cross-linked into 
one bundle, e.g. by the Ena/VASP proteins or the 'tip com- 
plex' 1 13], the aggregate can still be described as a single 
fiber obeying Eqs. 12151 up to a rescaling 5 — > ^, k 0Tli0 ff —> 
n k n,off and n —> n 2 k, which amounts to assuming that the 
fibers are strongly cross-linked. As a result of bundling both 
stalling and run-away occur at larger membrane forces. 

Fig. 4 shows how the dynamic 'phase diagram' of Fig. 
2 gets modified by bundling. Here we have taken n = 10, 
which may be a reasonable assumption for small filopodia 
in vivo |7|]. Solving Eqs. 12151 shows that C decreases with 
n, so that the stalling regime widens (Fig. 4). More im- 
portantly, a stiffer fiber, or a bundle with larger n, will de- 
form the membrane more. As a result if n is large enough, 
and #o is close to 90°, non-linear membrane deformation into 
tubular spicules 0] may occur before run-away. We adopt 
a rough test for when spicules will be formed that is related 
to the breakdown of the linear membrane Hamiltonian (1), 
specifically that the maximum average membrane gra dient 
exceeds unity {\V ±u\ max ) > 1. It can be shown 1 1 Qfl that 
this spiculation condition can be translated to a condition on 
the average membrane displacement, namely that it exceeds 

(zm) — d + JX dl ° g dd tot ^ 9 wn ich can be calculated from ^ 
(z m ) — d=^ ^/K m a. Inserting appropriate physical values 
we estimate that spicules form only if the (maximum) mem- 
brane displacement is at least ~ 100 nm or, equivalently, that 
the normal membrane force is ~ 80 pN M14I1 . Note that the 
boundary between spicule formation and fiber run-away is not 
sensitive to the polymerisation rate, unlike the boundary be- 
tween spicule formation and stalling. Thus, for physiological 
parameters, there is only a finite range of incidence angle for 
which a bundled actin fiber can form a filopodium. 

We know of at least two other mechanisms beyond fiber 
bundling which the cell may employ to avoid excessive run- 
away of actin fibers: (i) capping proteins Q], (ii) tethering of 
fiber tips to the membrane 

It is also interesting to consider a series of in vitro experi- 
ments aimed at controlling cell motility via protein expression 
fl^l in lights of our results. These experiments showed that 
lamellipodia only form when capping proteins are present in 
the cell, so that they can be thought of as a dense network 
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FIG. 4: Phase diagram for a bundle of n = 10 actin fibers, with pa- 
rameters obtained by scaling those used for Fig. 2 to this bundle size 
(see text). The insets are 3dMC simulations obtained with parameter 
values as in Fig 2 except the spicule for which y/AS = 0.1. 



of short and mostly artificially stalled fibers. Depleting the 
cell of capping proteins, without decreasing the amount of 
Ena/VASP proteins, results in massive growth of filopodia, 
akin to spicules. Reducing the number of Ena/VASP proteins 
as well, leads on the other hand to the formation of mem- 
brane ruffles, which may be associated with extensive fiber 
run-away, in agreement with what we predict for actin fibers 
polymerising with cell-like conditions. All the states in our 
phase diagram can thus be obtained in a cell by regulating the 
concentrations of these key proteins. 

The validity of our treatment is not limited to single and 
bundles actin fibers. Microtubules and sickle hemoglobin 
fibers may have parameters closer to 5 ~ 0.5 nm and n ~ 
0.1 — 1 mm. If we adopt d ~ 1 /im, comparable to the size 
of the cell, stalling occupies a smaller region in the dynamic 
'phase diagram' , and run-away and spicule formation a larger 
region. This is in agreement with experiments |5|]. 

Nonlinear membrane deformations, leading to spicule (or 
filopodia) formation cannot be described quantitatively us- 
ing Eqs. 12151 To check our approximate analytical threshold 
for the onset of spicule formation, we have performed full 3- 
dimensional Monte-Carlo dynamic simulations, in which both 
membrane and fiber are discretised into connected 'beads'. A 
time step in the dynamic evolution consists of attempting to 
move the position of all beads in both the fiber and the mem- 
brane via local moves (see Refs. ISlllal for details on the algo- 
rithm). A local move is accepted according to the Metropolis 
test with discretised versions of Eq.[T] 

Examples of run-away state and spicule, found via our sim- 
ulations, are shows as insets in Figs. 2 and 4. (Parameters 
are given in the legends.) For 0q = 90°, we confirm that 
there is a crossover between run-away and spiculation. For 
VAS ~ 0.2 this occurs when the persistence length of the 
polymers exceeds 1.5 — 2 /im, in good agreement with 1.25 
fim found using our Gaussian approximation. In few cases for 
lower l p values the fiber does not run away, but is 'trapped' by 
the membrane fluctuations giving rise to a spicule. We also 
find that the run-away state (see Figs. 2 and 4) gives rise to 
a bump in the membrane which may account for the observa- 
tion of ruffles 11511 . Finally, if we require that a fiber does not 
run-away for 10 8 Monte-Carlo steps, 3dMC simulations pre- 



dict that the crossover to the run-away state in Fig. 1 occurs 
for C ~ 85° for k on /k Q ^ = 100, which is in good agreement 
with C ~ 86° found via Eqs. EE 

In conclusion, we have solved the dynamic equations of 
motion for a growing semi-flexible polymer, incident onto a 
fluctuating soft membrane at an angle, and found a dynamic 
transition between fiber run-away and stalling. Given the ki- 
netic polymerisation rate of actin fibers in the cell, we predict 
that single actin fibers should bend and run away from the 
membrane, and would not contribute to cell motility. This is in 
agreement with in vitro investigations under unphysiological 
conditions in which the behaviour of cells lacking bundling 
and capping proteins was investigated. Moving cells seems to 
have taken 'precautions' to stop this phenomenon from hap- 
pening excessively, by bundling or capping fibers to render 
them stiffer. Our phase diagram can be generalised to the 
case of an actin bundle, and we show that in this case a third 
regime appears, in which the fiber deforms the membrane in a 
non-linear fashion and spicules appear, which are akin to the 
filopodial protrusions observed in motile cells. It is intrigu- 
ing that physiological values of the thermodynamic parameter 
log k on /k fi can be seen to lie in, or close to, the narrow range 
of values for which actin filipodia can explore all three differ- 
ent destinies. Finally, we predict that microtubules and sickle 
hemoglobin fibers should be more likely to form spicules or 
to run away than to stall, in agreement with observations. 
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